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Abstract 

A general model for the fragmentation of a two-component system (e.g. 

protons and neutrons) is proposed and solved exactly. The extension of this 

model to any number of components is also shown to be exactly solvable. A 

connection between this models and the permutation group is discussed. The 

notion of isotopic equivalence is defined in order to evaluate the equivalence of 

these models to earlier one-component models. All the one-component models 

considered in earlier papers are shown to be equivalent to a particular class of 

two-component models. A simplified model applicable to the case of nuclear 

fragmentation is introduced and analyzed. Modifications to this model to 

include effects such as pairing and Coulomb interactions are discussed. 
25.70.Pq, 05.30.Ch 
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I. INTRODUCTION 



The fragmentation of mesoscopic systems, particularly heavy nuclear systems, has at- 
tracted considerable interest over the past several years [|ll-p3|. A number of statistical 
models have been proposed to explain the phenomena. In these models, the final state of 
the fragmentation of a nucleus of A nucleons is represented by the vector n = [rii, . . . , ua), 
where rij is the number of fragments with j nucleons. The set of all such states or partitions 
of A is denoted by 11^, each of which must satisfy J2f=ijnj = A. 

In actuality, the final state of a nuclear fragmentation process is a more complicated 
object. Fragments are formed with varying numbers of neutrons and protons, which are 
distinguishable objects. For this reason, a representation of the final fragmentation state 
of such a two-component or isotopic system should minimally contain information on the 
neutron and proton content of each fragment. In this representation, the final state of 
the fragmentation of a nucleus containing Z protons and N neutrons is given by a vector 
of {Z + 1){N -|- 1) — 1 integers, n = (noi, . . . ,noN,nio, . . . ,niN, ■ ■ ■ ,nzN), where njk is the 
number of fragments with j protons and k neutrons. The set of all such states or partitions of 
A = Z + N is denoted by II^jv, each of which is subject to the two constraints Yl,jkj^jk = 
Hjk krijk = N. 

In general the fragmentation of an A^-component system oi A = YljLi particles where 
Aj is the number of indistinguishable particles of type j can be described by a vector 
of Ilf=M, + 1) - 1 inte gers, n = (noo...oi, • • • , '^Ai...Ajv), where nki,„kj^ is the number of 
fragments with kj particles of type j, j = 1, . . . ,N. The set of all such partitions is denoted 
by n^i..,Ajv5 sach of which is subject to the constraints J2ki...kM kj^ki...kN = J = 1; ■ ■ • ; ^• 

There are two principle reasons why one-component models are primarily used in the 
statistical modeling of nuclear fragmentation. The first reason is simplicity. One-component 
models are fundamentally simpler than two-component models, and the exact solution for a 
class of one-component models has been known for some time. In addition, as Sneppen 
noted, the handling of the full two-component system increases the number of allowable 
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final states dramatically (see table |). This makes it impractical to model these systems 
by standard Monte Carlo techniques. As these methods have been the primary technique 
used in evaluating realistic models which do not admit exact solutions, this discourages 
approaching the problem this way. 

The second reason is sufficiency. A one-component model can always model the physics 
if the underlying isotopic states are inaccessible to the observer. For example, typically 
only the charge distribution of the fragments is measured in a fragmentation experiment. 
The isotopic distribution for any particular charge Z is not recorded. So a one-component 
model distinguishing protons is sufficient to model the observables measured. The weight 
for a particular charge fragment distribution is then given by summing the two-component 
weights of all the states which have an equivalent configuration of protons. Since this 
weight is a complicated sum of two-component weights, there is no guarantee that the one- 
component weight will be simple. Indeed, this is usually not the case, as we will discover in 
section P^ . 

However, there are good reasons to prefer modeling nuclear systems with a two- 
component model, especially if such a model is amenable to analytic solution. For instance, 
any choice of statistical weight based on the fundamental statistics of the constituent par- 
ticles must by necessity be based on a two-component model. Also, unlike one-component 
models, definite predictions about the isotopic distribution of fragments could be made for a 
two-component model. Additionally, nuclear effects such as proton-neutron symmetry, pair- 
ing and Coulomb interactions could be modeled, effects which are difficult (if not impossible) 
to treat in a one-component model. 

This paper shows that a particular class of two-component models are in fact exactly 
solvable by recursion. Realistic models lie within this class, and their solution obviates 
the need for Monte Carlo and allows for the exact evaluation of all expectation values. A 
discussion of the relation of partitions to the permutation group Sa provides a closed form 
solution for a specific two-component model, as well as a proof that all the one-component 
models of the same form as those in ref. @] are equivalent to a particular two-component 
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model. 



II. ISOTOPIC MODELS AND THEIR SOLUTION 



Earlier papers discussed the one-component model where each state n G Hz was 
given the weight 

Wz{n,x) _J_Ax7 
Pz{n, x) = = —J— 11 ^ (1) 

Zz{x) Zz{x) fj-^ rijl 

The partition functions can be obtained by recursion, 

Zz{x) = WjxjZz-,{x) (2) 
^ j=i 

and all the ensemble averages are simple functions of the partition function. 

In analogy with the one-component model, we can define a canonical weight for n G Hzn 

by 



WzN(n,x) 1 -^x 



PzN{n,x) = = ^^11-^ (3) 



where x is a parameter vector which determines the fragmentation behavior. The l/njkl are 
the expected Gibbs factors which arise from the exchange statistics of identical fragments. 
The solution of this system can be derived from the following property of the partition 
function: 

dZzN ry (A\ 

-Q^ = ^^-^^-^•(^) (4) 

where Zjk{x) = when j < or /c < 0. 

From this property any moment can be calculated in terms of the partition functions. 
For example, 

^ Zzn{x) dx,k ' Zzn{x) ^ ^ 

or in general. 
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n^Pjk ^Z-J2jPjk,N-J2 kpjk (^) 
jk Zzn[X) 

where p is a vector of integers, p — J2jkPjk ^-nd the faUing factorial [z\k is given by [z]k — 
{z - k + l)[z]k-i, with [z]o = 1. 

Of course, the partition functions are needed to compute these expressions. The con- 
straints J2jk j'l^jk — Z, J2jk krijk = N hold for all configurations, so they must also hold 
for the expectation values. This allows us to express the partition function as a recursively 
defined function: 

Zzn{x) = ^Yl^XjkZz-j,N-k{x) 
jk 

Zzn{x) = J^Y1 kxjkZz-j,N-k{x) . (7) 

jk 

Since Zoo{x) — 1, these two relations can be used to iteratively construct the set of partition 
functions for any parameter vector x. 

It should be noted that the technique used to solve this model applies equally well to 
any number of components. The A'^-component system of ^ = ^f=i particles with the 
weight 

Pm...aMx) = - n (8) 

^Ai...AN\X) ki...kN ^ki...kN- 

has a recursively defined set of partition functions determined by the N equations 

Zai...aA^) = ^ II kjXki...kNZA,-ki,...,Aj,-kM{^) (9) 

•J fci...fejv 

for j = 1, . . . , where ^oo...oo(a^) = 1- AH the ensemble averages can be determined solely 
from these partition functions and the parameter vector x. 



III. RELATION TO THE PERMUTATION GROUP 



It is useful to consider an alternative representation of these models. Consider the 
set of all permutations on A objects, the symmetric group 5"^. Typically a permutation 
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p G Sa is represented by a vector of integers p = {pi, ■ ■ ■ ,Pa) where pj is the result of the 
permutation acting on j, i.e. j ^ pj. We can also represent p by its cycle decomposition. 
In this representation, the permutation is a set of disjoint cycles of varying length. Each 
cycle contains all the numbers which permute among themselves under the action of the 
permutation. 

If we identify each element in the permutation vector with a particle, we can consider 
a cycle in the cycle decomposition of a permutation as representing a cluster of particles. 
This defines a map from the set of permutations Sa onto the set of partitions 11^ by cycle 
decomposition. For example, table || details this for A = A. The number of permutations 
which map to a particular partition of A was computed by Cauchy, and is given by 

Since this is of the same form as the weight given by the one-component partition model, 
we see that the model given by the permutation weight 



is equivalent to the partition model given by eq. (|^) with xj given by 

X. = ^ (12) 

J 

We can also map the set of permutations Sa onto the set of partitions II^at where 
A = Z + N in an analogous way. Identify elements 1, . . . , Z in the permutation as protons, 
and elements Z + 1, . . . , Z + N as neutrons. Each cycle in the permutation is mapped to 
a cluster by counting the number of "neutrons" and "protons" in that particular cycle. As 
an example, this is done for Z = 2, = 2 in table The number of permutations which 
map in this way to a particular n G 11^^^ is given by the two-component generalization of 
Cauchy' s formula 

M.in:>^z^.N^.n^(^l±±piy' (13, 



which satisfies the sum rule J^hgUzn = {Z + N)\. Since this number is of the same 

form as eq. (H), the permutation model given by 

PzNip,^ = ^^Y[^T^'^ (14) 

Zzn{z) 

is equivalent to the model given by eq. (0) with 

+ .... 

This can be generalized to any number of components. The iV-component generalization 
of Cauchy's number for A = J2f=i is given by 

M.(ri) = n -4,,! n %rVr 

j=l k^...kj, ^k,...kN- \ Ki\---kN\ J 

and the permutation model 

pm...aAz) = ^ ' n (17) 

is equivalent to the partition model given by eq. (H) with Xki...kj^ given by 

Using the two-component model as defined over the permutation group allows us to 
construct a closed-form solution for the partition function of a particular model. Consider 
eq. (|14D with Zjk = x, where a; is a parameter characterizing the degree of fragmentation 
of the system. Each permutation contributes x"^ to the partition function, and there are 
\S^^\ permutations with m fragments, where is the Stirling number of the first kind. 
Therefore Zzn{z) = x{x + 1) . . . {x + A — 1). This implies that the partition model with 

U + k~ 1)! X (j + k\ 

has its partition function given by the closed-form expression 

1 T(x + Z + N) 



The expected number of fragments with j protons and k neutrons is 

(n,,) = X (^^^ (J^^ B{j + k,x + A-j-k) (21) 

where B{x,y) is the Euler Beta function, B{x,y) = T{x)T{y) /T{x + y). The expected 
number of fragments of charge j, Uj = Y.k=o''^jk is given by 

{n,)=x(^^'jB{j,x + Z-j) (22) 

which can be shown by applying the Norlund formula [x + y]" = J2k=o {k}[^]''[yT~^ 
eq. (pTD, where the rising factorial is given by [x]^ = {x + k — and = 1. 

Notice that eq. (p^) is identical to the one-component expectation value for the model 
given by eq. (|l]) with xj = x/j (see 0). This is not a coincidence. In fact, the relation of 
the partition models to the permutation group allows us to prove in the next section that 
there exist two-component models of the form of eq. (|^) for every one-component model of 
the form of eq. (|T|) such that the results of the one-component model are preserved when 
a summation over isotopes is performed. In this case, the two-component model given by 
eq. ([ToD is equivalent to the one-component model given by xj = x/j. 

IV. ISOTOPIC AND PERMUTATION EQUIVALENCE 



As mentioned in section |T|, two-component models are sometimes equivalent to one- 
component models. In this section, we make this notion precise as well as show under what 
conditions this is true. 

A two-component model is said to be isotopically equivalent to a one-component model 
if for every function /(n), the one-component ensemble average {f)z = J2n f{^)Pz{n) is 
equal to the two-component ensemble average {f)zN = J2n' f{^{^'))f"zNi'^') where nj{n') = 
Y^k=o^'jk- example, {nj)z = Y^k=o{'^3k) zn must hold if two models are considered to be 
isotopically equivalent. As mentioned in the introduction, this is true if and only if 

Pz{n)= E P'zNin') (23) 

n'&UzNir!,) 



where YlzNin) denotes all the partition vectors ft' e Hzn such that rij — J2k=o''^'jk- 

If we consider deriving the models over the permutation group, we can define the notion 
of permutation equivalence of P'z'N'ip') to Pzn{p) by the condition 

Pzn{p)= E Pz'N'ip') (24) 

where Z' > Z , N' > N and Sz'+n'{p) is the set of all permutations p' e Sz'+n' which give 
p by the following construction. Decompose p' into its cycle representation. Eliminate the 
numbers Z + 1, . . . , Z' and Z' -\- N -\- 1^ . . . ^ Z' -\- N' from the cycles. Renumber the elements 
Z' + 1, . . . , Z' + as Z + 1, . . . , Z + to recover the permutation in Sz+n- In other 
words, if k is an element to eliminate, find j such that pj — k, and modify the permutation 
vector so that pj — pk- This removes the element from the permutation. Do this for all 
the numbers to eliminate, resulting in a new permutation isomorphic to a permutation in 
Sz+N- For example, p' — (1, 5)(3)(4, 6, 2) e becomes (1, 5)(4, 2) by eliminating 3, 6. After 
renumbering this becomes p — (1, 4) (2, 3) e 5*4. 

Permutation equivalence of Pzn{p') to Pzoip) implies isotopic equivalence of 
M2{n'{p'))Pzj^{n'{p')) to M2{n{p))Pzo{'n{p)), assuming the permutation weights are only 
functions of n{p), i.e. all permutations which have the same partition configuration have the 
same weight. One particular permutation weight which satisfies this condition is the unnor- 
mahzed weight Wzn{p, z) — Hj-jt z]^t^^ ■ ^^'^ show that this weight with A — Z + N 

is permutationally equivalent to the Z,N + 1 weight with the same Zjk if Zjk — Zj. This 
implies by induction that the ZO model is permutationally equivalent to the ZN model, and 
therefore that the two weights corrected by the Cauchy factor are isotopically equivalent. 

For any permutation p & Sa, there are A + 1 permutations in SA+iip)- One of these per- 
mutations has the element A+1 in its own cycle. It has weight Wz,N+iip', z) = Wzn{p, z)zoi. 
A of the permutations are hybrids, with the extra element combined with a cycle from p. 
Suppose the element is in a cycle of j neutrons, k protons. The weight of that cycle is 
Wz,N+i{p' , z) = Wzn{P: z)zj^k+i/ Zjk, since the introduction of the additional element has 
reduced rijk by one, but increased rij^k+i by one. There are j + k possible places for the 

9 



element to be inserted into that cycle, and there are njk cycles of that type it can combine 
with. So the following is true 

^ Wz,N+l{p') . ^ . .^.x 

where fjk = (j + k)zj^k+i/zjk- If -Zjfc is independent of k, this reduces to 2:01 + + Z. In 
this case, we can multiply both sides of the equation by Wzn{p) and sum over Wp G Sa to 
arrive at 

Zz,N+ii^ = Zzn{z){zoi + Z + N). (26) 

This implies eq. (|2^ ) when substituted back into eq. (^). So under these conditions, Pzn 
is equivalent to Pz,n+i- Therefore, Pzo is equivalent to Pzn by induction and transitivity, 
and eq. (P is equivalent to eq. @ for Xjk given by 

Xjk = Xj ^ ^ (27) 

The partition function for this two-component model is 

ZMi) = (28) 

Notice that at no point was it asserted that these are the only two-component models 
that are isotopically equivalent to the one-component models. In the previous discussion, 
the key property was the independence of J2jk fjknjk from n. Perhaps we could satisfy 
this with choices other than fjk = j + k. In fact, it can be proven that this condition is 
satisfied only by fjk = aAj/Z + (1 — a)Ak/N where a is arbitrary. However, a G [0, 1] 
to guarantee nonnegative probabilities in the statistical model. Given ZjQ, we could then 
construct permutationally equivalent two-component models specified by a and the recursion 
relation 

A ( j k\ 

and evaluate the equivalent partition weights using eqs. (H), ([151) . 
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If we require that the models be isotopically equivalent for any choice of N, these are 
the only models given by eqs. (|I]), (|^) that are isotopically equivalent. This is true because 
isotopic equivalence for any implies permutation equivalence of to + 1. Permutation 
equivalence requires J2jk fjknjk to be independent of n, which can be shown to be true only 
if Zjk is of the form given above. Notice that this condition is very similar to Kingman's 



non-interference condition on one-component partition weights p4[. We require that the 
neutrons do not "interfere" with the overall proton distribution, just as in one-component 
models Kingman required that the addition of new objects does not interfere with the 
probability structure of the original set of objects. 

This technique can be generalized to construct an A^ + l-component model that is equiva- 
lent to an A^-component model. For example, if the A^-component weight is given by eq. (^, 
then an equivalent A^ + 1-component model is given by the parameter vector defined by 

' (30) 

These results provide a simple way of constructing two-component models from the one- 
component models explored in earlier papers. For example, one could apply eqs. (P7D, (^S|) 
to the chain model, Xj = x, to construct a two-component extension of the chain model as 



developed by Gross, et. al. 10 



V. APPLICATION TO NUCLEAR FRAGMENTATION 

The two-component analogs of one-component models given by eq. (^) could be taken 
as the starting point for a fully isotopic model of fragmentation. These models would contain 
no surprises in their one-component expectation values, and would provide predictions for 
the isotopic behavior. 

This is perhaps too strong a requirement. Equation (|29| ) constrains the set of possible 
partition weights, forbidding a number of physically interesting choices for xjk- For example, 
if the symmetric term exp{—as{j — kY/ {j + k)} appeared in the parameter zjk, the resulting 
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parameter vector could never satisfy eq. (pQl). Since such a symmetry term should appear 
due to the symmetry term in the binding energy of nuclei, we should not expect realistic 
models to be isotopically equivalent to the models considered in earlier papers. 

With that in mind, we start with a simplified model of nuclear fragmentation specified 
by the following choice for Xjk 

X (3 + k\ j asjj - kf \ , . 

"^'=(7wl 7 j^'^n-^T^r '''' 

We will find it appropriate sometimes to modify this when j + k = 1, i.e. to avoid the 
exponential suppression factor e^"" for monomers. Let us specialize to the case r = 1 and 
consider the limiting cases. When a<( — > 0, the model becomes the two-component analog of 



the model discussed in section |ITT| . The partition function has a closed-form solution, which 
is given by eq. (]20|). When oo, Xjk is proportional to a Kronecker delta function and 

the model reduces to the one-component model specified by parameters 

where k indexes both proton and neutron number, since only symmetric fragments with 
j = k are allowed. In this case, the partition function also has a closed-form solution 

X r{x + 2A) 
= A^.Tix + A + l) 

As an aside, note that all models with Xk = ^{^^^ for n = 1,2,3,... have closed-form 
solutions given by Za{x) = {x/A\){r{x + nA)/r{x + (n — 1)^4 + 1)). Models with interme- 
diate values for are not as simple as these two limits, but can still be computed easily. 
Figures |l|, ^ show the prediction of this model at x = 1 and various as for a symmetric case 
Z = N = 50 a.s well as an asymmetric case Z = 50, = 75. These figures display the total 
charge distribution {nj) = I]fcLo("-jfc)) &s well as the distribution of isotopes for a particular 
charge, specifically carbon, i.e. {njk) for j = 6. 

From the total charge distribution figures, we see that increasing increases the num- 
ber of large fragments. This is due to the exponential suppression of asymmetric fragments, 
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which make up a large fraction of the partition space, especially for smaller fragments. The 
isotopic distribution figures reveal that the distribution of isotopes is essentially Gaussian. 
As seen earlier 0, an initial asymmetry in the relative numbers of protons and neutrons 
yields a similar asymmetry in the distribution of fragments. This asymmetry can be dimin- 
ished by increasing and removing the exponential suppression from monomer parameters 
(xoi, xio), as was done in the figures. Indeed, when as — > oo, such a model would allow only 
symmetric fragments and monomers. 

To this simple model, we can add other effects, eventually including all the terms in the 
semi-empirical mass-formula. For example, nucleon pairing interactions can be added using 
the term exp{apPjk/ {j + k)} where pjk is for even-odd nuclei, —1 for odd-odd nuclei, and 
+1 for even-even nuclei. This can be done for all intrafragment interactions. Interfragment 
interactions do not fit the form of eq. (j^) and must be handled in some kind of mean- 
field manner in this approach. Coulomb interactions are the only persistent long-range 
interactions in the thermalized phase. Suppose each fragment interacts with all the others 
via the Coulomb interaction. This contributes an interaction energy given by 

El = ^Yl ^Jkinim - 5jl5km)E]klm (34) 
jklm 

If we choose -E/fc^m = 3^^^ I ('^) ^^e Coulomb interaction, then 

Ei = ^AZ'-Y.j\k) (35) 

which suggests we include interfragment Coulomb interactions by appending the term 
exp(acj^) to Xjk- We can estimate ac by assuming (r) ^ V^^^ ~ tqA^^^. 

We postpone investigating such a realistic choice for Xjk for a later paper. 

VI. CONCLUSION AND SUMMARY 

This paper established a non-trivial exactly solvable isotopic model of fragmentation. 
The model is specificed by a set of parameters xjk, where Xjk includes all the physical 
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parameters pertinent to the fragmenting behavior of clusters of j protons, k neutrons. This 
model is related to the permutation group and as a result, it was proven that for certain 
choices of Xjk it directly reduces to a one-component model when a summation over isotopes 
is performed. The usefulness of this model in nuclear fragmentation is manifest, and a 
simplified model analyzed here contains most of the salient features of the phenomenon. A 
more realistic model is also discussed. 
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FIGURES 

FIG. 1. Overall (a) and Isotopic (b) distribution of fragments for Z = 50, AT = 50, x = 1, r = 1 

and various Ug. Isotopic distribution is for carbon, Z = 6. 

FIG. 2. Overall (a) and Isotopic (b) distribution of fragments for Z = 50, = 75, x = 1, r = 1 
and various a^. Isotopic distribution is for carbon, Z = Q. 
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TABLES 



TABLE I. Number of partitions for Z protons and N neutrons. 

Z N Pz Pn+z Pzn 

10 10 42 627 59521 

20 20 627 37338 3026 73029 

30 35 5604 20 12558 231 50477 43167 

40 50 37338 566 34173 6 00944 62493 61633 

50 70 9 66467 18443 49560 26953 04447 12166 27712 
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TABLE II. The twenty-four permutations of S4 and their partition representations. The first 
two columns give the typical representation and cycle decomposition of the permutations. The third 
column gives the equivalent one-component partition vectors n = (n-i, n2, ns, 714). The last column 
gives the equivalent two-component partition vectors n = {noi,no2,nio,nii,ni2,n2o,n2i,n22) for 

Z = 2. N = 2. 



p 


P 


n4 


n22 


(1,2,3,4) 


(1)(2)(3)(4) 


(4,0,0,0) 


(2,0,2,0,0,0,0,0 


(1,2,4,3) 


(1)(2)(3,4) 


(2,1,0,0) 


(0,0,2,0,0,1,0,0 


(1,3,2,4) 


(1)(2,3)(4) 


(2,1,0,0) 


(1,0,1,1,0,0,0,0 


(1,3,4,2) 


(1)(2,3,4) 


(1,0,1,0) 


(0,0,1,0,1,0,0,0 


(1,4,2,3) 


(1)(2,4,3) 


(1,0,1,0) 


(0,0,1,0,1,0,0,0 


(1,4,3,2) 


(1)(2,4)(3) 


(2,1,0,0) 


(1,0,1,1,0,0,0,0 


(2,1,3,4) 


(1,2)(3)(4) 


(2,1,0,0) 


(2,1,0,0,0,0,0,0 


(2,1,4,3) 


(1,2)(3,4) 


(0,2,0,0) 


(0,1,0,0,0,1,0,0 


(2,3,1,4) 


(1,2,3)(4) 


(1,0,1,0) 


(1,0,0,0,0,0,1,0 


(2,3,4,1) 


(1,2,3,4) 


(0,0,0,1) 


(0,0,0,0,0,0,0,1 


(2,4,1,3) 


(1,2,4,3) 


(0,0,0,1) 


(0,0,0,0,0,0,0,1 


(2,4,3,1) 


(1,2,4)(3) 


(1,0,1,0) 


(1,0,0,0,0,0,1,0 


(3,1,2,4) 


(1>3,2)(4) 


(1,0,1,0) 


(1,0,0,0,0,0,1,0 


(3,1,4,2) 


(1,3,4,2) 


(0,0,0,1) 


(0,0,0,0,0,0,0,1 


(3,2,1,4) 


(1,3)(2)(4) 


(2,1,0,0) 


(1,0,1,1,0,0,0,0 


(3,2,4,1) 


(1>3,4)(2) 


(0,0,0,1) 


(0,0,0,0,0,0,0,1 


(3,4,1,2) 


(1,3)(2,4) 


(0,2,0,0) 


(0,0,0,2,0,0,0,0 


(3,4,2,1) 


(1,3,2,4) 


(0,0,0,1) 


(0,0,0,0,0,0,0,1 


(4,1,2,3) 


(1,4,3,2) 


(0,0,0,1) 


(0,0,0,0,0,0,0,1 


(4,1,3,2) 


(1,4,2)(3) 


(1,0,1,0) 


(1,0,0,0,0,0,1,0 


(4,2,1,3) 


(1,4,3)(2) 


(1,0,1,0) 


(0,0,1,0,1,0,0,0 
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(4.2.3.1) (1,4)(2)(3) (2,1,0,0) (1,0,1,1,0,0,0,0) 

(4.3.1.2) (1,4,2,3) (0,0,0,1) (0,0,0,0,0,0,0,1) 
(4,3,2,1) (1,4)(2,3) (0,2,0,0) (0,0,0,2,0,0,0,0) 
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